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ABSTRACT 

The unsplittable flow problem is one of the most extensively 
studied optimization problems in the field of networking. An 
instance of it consists of an edge capacitated graph and a 
set of connection requests, each of which is associated with 
source and target vertices, a demand, and a value. The 
objective is to route a maximum value subset of requests 
subject to the edge capacities. It is a well known fact that 
as the capacities of the edges are larger with respect to the 
maximal demand among the requests, the problem can be 
approximated better. In particular, it is known that for 
sufficiently large capacities, the integrality gap of the corre- 
sponding integer linear program becomes 1 + e, which can 
be matched by an algorithm that utilizes the randomized 
rounding technique. 

In this paper, we focus our attention on the large capac- 
ities unsplittable flow problem in a game theoretic setting. 
In this setting, there are selfish agents, which control some 
of the requests characteristics, and may be dishonest about 
them. It is worth noting that in game theoretic settings 
many standard techniques, such as randomized rounding, 
violate certain monotonicity properties, which are impera- 
tive for truthfulness, and therefore cannot be employed. In 
light of this state of affairs, we design a monotone determin- 
istic algorithm, which is based on a primal-dual machinery, 
which attains an approximation ratio of ^§-j- , up to a dispar- 
ity of e away. This implies an improvement on the current 
best truthful mechanism, as well as an improvement on the 
current best combinatorial algorithm for the problem under 
consideration. Surprisingly, we demonstrate that any algo- 
rithm in the family of reasonable iterative path minimiz- 
ing algorithms, cannot yield a better approximation ratio. 
Consequently, it follows that in order to achieve a mono- 
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tone PTAS, if exists, one would have to exert different tech- 
niques. We also consider the large capacities single-minded 
multi-unit combinatorial auction problem. This problem is 
closely related to the unsplittable flow problem since one can 
formulate it as a special case of the integer linear program 
of the unsplittable flow problem. Accordingly, we obtain a 
comparable performance guarantee by refining the algorithm 
suggested for the unsplittable flow problem. 

Categories and Subject Descriptors 

F.2 [Theory of Computation]: Analysis of Algorithms 
and Problem Complexity 

General Terms 

Algorithms, Economics, Theory 

Keywords 

Mechanism design, approximation algorithms, combinato- 
rial and multi-unit auctions, primal-dual method 

1. INTRODUCTION 

The problems. We study the unsplittable flow problem. 
As input to this problem, we are given a directed or undi- 
rected graph G — (V,E), such that n — \V\, m — \E\, and 
every edge e £ E has a positive capacity c e . An additional 
ingredient of the input is a set 1Z of connection requests, in 
which every request r G 1Z is characterized by a quadruple 
(s r ,t r ,d r ,v r ) such that s r and t r are the respective source 
and target vertices of the request, d r is the positive demand 
associated with the request, and v r is the positive value or 
profit gained as a result of allocating the request. The objec- 
tive is to select a maximum value subset of requests S C 1Z, 
along with a path for each selected request, so that all the 
requests in S can simultaneously route their demand along 
the corresponding path, while preserving the capacity con- 
straints. Denoting by B — min e {c e }/max r {d r } the ratio 
between the minimal capacity of an edge and the nraxinral 
demand among the requests, the problem is referred to as 
the B-bounded unsplittable flow problem. Since one can nor- 
malize both the demands of the requests and the capacities 
of the edges, the B-bounded unsplittable flow problem can 
be equivalently defined to have d r € (0, 1] for every request 
r, and B = min e {c e }. Note that we shall use the latter 
definition throughout this paper. 

We also consider the single-minded multi-unit combina- 
torial auction problem. This problem is closely related to 



the unsplittable flow problem since one can formulate it as 
a special case of the integer linear program of the unsplit- 
table flow problem. An instance of it consists of a set U of 
m non-identical items, such that item u £ U has a positive 
integer multiplicity c u . The input also consists of a set TL of 
requests, in which every request r £ 1Z is characterized by a 
pair (U r ,v r ) such that U r C U is an items bundle, which is 
the demand associated with the request, and v r is the pos- 
itive value gained as a result of allocating the bundle. The 
goal is to select a maximum value subset S C 1Z, so that 
every item u £ U appears in at most c u bundles of requests 
in S. Denoting by B = min u {c u } the minimum multiplic- 
ity of an item, the problem is referred to as the B-bounded 
multi-unit combinatorial auction problem. 

The setting. In the present paper, we study the f2(lnm)- 
bounded unsplittable flow problem, and the Q(lnm) -bounded 
multi-unit combinatorial auction problem from a mechanism 
design [15] point of view. In this game-theoretic setting, 
some characteristics of the requests, which are henceforth 
referred to as the type of the requests, are controlled by self- 
ish agents. An agent is selfish in a sense that it might declare 
a fallacious type in order to manipulate the algorithm in a 
way that will maximize its own utility. Our goal is to design 
mechanisms, which are referred to as incentive compatible or 
truthful, in which each agent's best strategy is always to re- 
veal the true type of the request that it controls, regardless 
of the other requests types, and regardless of the way that 
the other agents decide to declare their requests types. In 
particular, we aim to devise monotone algorithms which are, 
roughly speaking, equivalent to truthful mechanisms. Note 
that in the unsplittable flow problem, the type of a request 
is its demand and value, whereas in the multi-unit combi- 
natorial auction problem the type of a request is its value. 
Also note that the other characteristics of the request, e.g. 
the source and target vertices in the unsplittable flow prob- 
lem, are assumed to be known and thus, the agent cannot 
be untruthful about them. 

The motivation. One of the closely related problems to 
the unsplittable flow problem is the multicommodity flow 
problem. Since the multicommodity flow problem can be 
modeled by the relaxation of the integer linear program of 
the unsplittable flow problem, it may be considered as its 
fractional version. It is well known that the integrality gap of 
the integer linear program of the unsplittable flow problem 
becomes 1 + e when the ratio between the minimal capacity 
of an edge and the maximal demand among the requests 
is sufficiently large. Consequently, it is conjectured that in 
such case, the performance of algorithms for the fractional 
and integral versions would be similar. Specifically, since the 
multicommodity flow problem admits a monotone PTAS by 
combinatorial primal-dual based algorithms [9, 8], one may 
expect that an integral version of these monotone PTAS 
would yield a monotone PTAS for the unsplittable flow prob- 
lem. In the following, we refute this perception. In particu- 
lar, we design a primal-dual based monotone algorithm for 
the fi(lnm)-bounded unsplittable flow problem that attains 
the best possible approximation ratio with respect to any 
reasonable iterative path minimizing algorithm 1 , and is not 
a PTAS. Nevertheless, This algorithm still improves upon 
the best previous known result [7]. 



1.1 Our results 

Monotone deterministic algorithms. We describe a 
monotone deterministic algorithm based on a primal-dual 
approach for the fi(lnm)-bounded unsplittable flow prob- 
lem, which obtains an approximation ratio that approaches 
— 3j ~ 1.58. This result implies a corresponding incen- 
tive compatible mechanism. In addition, we show that the 
aforesaid algorithm can be attuned for the f2(lnm)-bounded 
multi-unit combinatorial auction problem, and attain a com- 
parable approximation ratio, i.e. — §j- approximation. These 
results improve over the approximation guarantee suggested 
by Briest et al. [7] for both problems, which approaches 
e»2.71. 

Deterministic lower bounds. We prove that any algo- 
rithm for the £l(lnra)-bounded unsplittable flow problem, 
which is part of the reasonable iterative path minimizing 
algorithms family, cannot yield an approximation guaran- 
tee that is better than ^§j — o(l). This implies, on the 
one hand, that the analysis of our algorithm is tight, and 
on the other hand, that to achieve a monotone determin- 
istic PTAS, if exists, one would have to employ different 
techniques. Additionally, we reinforce this inapproximabil- 
ity result by demonstrating that even if we ease the problem 
setting, e.g. assume that the minimal capacity of an edge is 
arbitrarily large, still no reasonable iterative path minimiz- 
ing algorithm can attain PTAS. Correspondingly, we also 
establish a lower bound of % on the approximation ratio of 
any reasonable iterative bundle minimizing algorithm for the 
n(lnm)-bounded multi-unit combinatorial auction problem. 

A deterministic (1 + e)-approximation algorithm. 

We study the f2(ln m) -bounded unsplittable flow with repe- 
titions problem, which is a variant of the fi(lnm)-bounded 
unsplittable flow problem in which one is allowed to satisfy 
every request multiple times using possibly multiple paths. 
In contrast with our prior findings, we demonstrate that this 
version admits a deterministic primal-dual based algorithm, 
which yields an (1 + e)-approximation. 

1.2 Related work 

The unsplittable flow problem, and its fractional variant, 
the multicommodity flow problem, has been given an exten- 
sive attention in recent years, both from an algorithmic point 
of view and from a game-theoretic one. These fundamental 
optimization problems model a diverse collection of appli- 
cations in network routing, parallel computing, and even in 
VLSI layout. Obviously, the fractional problem is easier. In 
particular, it is known to admit an optimal solution by lin- 
ear programming, and an (1 + e)-approximate solution by 
combinatorial algorithms, see e.g. [9, 8]. In contrast with 
the fractional problem, approximating the integral problem 
is hard. Guruswami and Talwar [11] have recently showed 
that the directed version of the problem is n ^'^-hard 
to approximate unless NP C BPTIME(n° (loglog " ) ), where 
B = 0(log nj log log n) is the minimal capacity of an edge. 
This result extended the prominent result of Guruswami 
et al. [10], which states that if B — 1, it is NP-hard to 
approximate the problem to within a factor of 0(n' 2 ~ € ). 
Respectively, when the graph is undirected, Andrews et 
al. [1] established an (logn) ^'^-hardness for any B — 
0(loglogn/logloglogn), under the assumption that NP $£ 



1 The family of reasonable iterative path minimizing algo- 
rithms is formally defined in Subsection 3.3. 



ZPTIME(- 



polylog(n) 



). Nevertheless, when B is sufficiently 



large, e.g. B — fi(lnm), the integrality gap of the inte- 



ger linear program of the problem becomes 1 + e, which 
can be matched by an algorithm that utilizes the random- 
ized rounding technique [17, 16, 18]. Unfortunately, this 
standard technique violates certain monotonicity properties, 
which are imperative for truthfulness and thus, cannot be 
directly used in the presence of selfish agents to obtain a 
truthful mechanism. Accordingly, until recently, the known 
truthful results for the f2(lnm)-bounded unsplittable flow 
problem only guaranteed approximation ratios that were at 
least logarithmic in the size of the graph [5, 6, 4]. Briest 
et al. [7] seem to have been the first to propose a constant 
factor approximation algorithm. Essentially, they designed 
a monotone primal-dual based algorithm, which was moti- 
vated by the novel work of Garg and Konemann [9] on the 
fractional problem, that achieves an approximation guaran- 
tee that approaches e. 

The research of the single-minded multi-unit combinato- 
rial auction problem, which is closely related to the unsplit- 
table flow problem, yielded similar results. Bartal et al. [6] 
showed that approximating the problem to within a factor of 
0(m 1 '( B+1 ' ) is NP-hard, where B is the minimum multiplic- 
ity of an item. Yet, when B = f2(ln 771), the integrality gap of 
the corresponding integer linear program becomes 1 + e. Ac- 
cordingly, Archer et al. [2], and Lavi and Swamy [12] devised 
truthful (1 + e)-approxfmation mechanisms. However, these 
mechanisms are truthful only in a probabilistic sense and 
hence, the best known deterministic truthful result for the 
fi(lnm)-bounded multi-unit combinatorial auction problem 
is by Briest et al. [7], which attains e-approximation. 

2. PRELIMINARIES 

In what follows, we present the notions of monotonicity 
and exactness, and then turn to describe a characterization 
that reduces the goal of designing truthful mechanisms to 
that of designing monotone and exact algorithms. Remark 
that the illustrated terms are presented in the context of the 
problems under considerations and hence, the keen reader 
may refer to Lehmann et al. [13] or Briest et al. [7] for more 
formal and comprehensive definitions of the underlying con- 
cepts. 

Definition 2.1. An algorithm A for the unsplittable flow 
problem is said to be monotone w.r.t. the demand and value 
of a request r € TZ, if it satisfies the following property: if 
algorithm A selects r when its demand is d r and its value is 
v r then algorithm A would have selected r if its demand was 
d r < d r , its value was v r >v r , and the demands and values 
of all the other requests were fixed. 

Definition 2.2. An algorithm A for the unsplittable flow 
problem is called exact, if it allocates the exact demand of 
every request selected, and does not allocate anything other- 



Theorem 2.3. ([13, 7]) If algorithm A for the unsplit- 
table flow problem is monotone and exact w.r.t. the demand 
and value of every request then there exists a correspond- 
ing truthful mechanism. In addition, this mechanism can be 
efficiently computed using algorithm A. 

Note that similar definitions can analogously be made for 
the single-minded multi-unit combinatorial auction prob- 
lem. The only exception is that the monotonicity property, 



and the characterization theorem are only defined with re- 
spect to the value of every request, i.e. the demand terms 
need to be cast off. 

3. UNSPLITTABLE FLOW PROBLEM 
3.1 The algorithm 

In this subsection, we devise a deterministic monotone al- 
gorithm for the f2(lnm)-bounded unsplittable flow problem, 
which achieves an approximation ratio that approaches — §j . 
Our algorithm is based on a primal-dual machinery. Accord- 
ingly, we present in Figure 1, the primal-dual formulation of 
the unsplittable flow problem. This will be later used to 
motivate the algorithm. 
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Figure 1: The integer linear program of the unsplit- 
table flow problem (top), and the dual of its relax- 
ation (bottom). Note that SV denotes the set of all 
the simple paths between s r and t r in G, S — {J re7i S r , 
and d s and v a denote the respective demand and 
value of path s, i.e. if s £ S r then d 3 — d r and v a = v r . 

Algorithm Bounded-UFP, formally described below, is a 
primal-dual based algorithm for the fi(lnm)-bounded un- 
splittable flow problem. Informally, the algorithm maintains 
the variables of the primal and dual programs, and in each it- 
eration selects to satisfy a request, which corresponds to the 
"most violated" constraint of the dual linear program. Fa- 
vorably, this reduces to finding a (normalized) shortest path 
in the graph G, whose edge weights correspond to the set of 
dual variables y e . It is worth noting that the algorithm, and 
part of its analysis is in the spirit of the algorithm suggested 
by Briest et al. [7] . 

We would like to note that since the path related variables, 
and the request related variables, i.e. the x 3 and z r variables 
respectively, play no role in the execution of the algorithm, 
lines 2, 3, and 12 are not regarded part of the algorithm. 
Nevertheless, we decided not to neglect them from the algo- 
rithm's description since they ease the analysis presentation. 



Algorithm 1 Bounded-UFP(e) 



Input: An accuracy parameter e G (0, 1] 
Output: A (request, path) pairs set W, which holds the 
requests to be allocated 

Let C be a list of all the requests, and let VV be an empty 

set 

for all r G C do z r = end for 

for all s G S do x s = end for 

for all e £ E do y e = — end for 

while (£ ^ and £ eSB c e J/ e < e 6 ^" 1 ') do 
for all r££ do 

Let p r be the shortest path between s r and i r in 
G with respect to the weights y e , and 

let |p r | = J2ee Pr Ve be its l en g tn 
end for 

Let f be the request, which minimizes — \p r \ with 
respect to every r G C 
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for all e G p? do y e = j/ e ■ e tBd f/ c ^ end for 
Add (f,pf.) to W, and remove r from £ 
Let a; p - = 1 and zt — Wf 
end while 

return W 



3.2 Analysis 

In this subsection, we will prove the following theorem. 

Theorem 3.1. For any e G (0,1], algorithm Bounded- 
UFP(jr) returns a feasible ((i+e) j§j)- approximate solution 



for thefl{ 1 - 



-bounded unsplittable flow problem, runs in 



polynomial-time, and is monotone and exact w.r.t. the de- 
mand and value of every request. 

Corollary 3.2. For all values e G (0,1], there exists 
a polynomial-time truthful ((i+<0— ^)- approximation mech- 
anism for the Q.(\nin) -bounded unsplittable flow problem, 
where every request's demand and value is unknown. 

We begin by introducing a notation that ease the analysis 
presentation: 

• Let x\, yl, and z r r be the respective values of the vari- 
ables x s , y e , and z r at the end of the i-th iteration of 
the algorithm, where i > 0. Mind that we regard the 
end of iteration as the beginning of the algorithm. 
Additionally, we let (y l , z l ) denote the set of dual vari- 
ables at the end of the i-th iteration. 



Let P(i) = Er 



(E s x\) be the value of the 



primal solution at the end of the i-th iteration, and let 
P be the value of the primal solution when the algo- 
rithm terminates. Notice that P is the sum of values 
of requests selected to be allocated by the algorithm, 
i.e. the outcome of the algorithm. In addition, we let 
AP(i) = P(i) - P(i - 1) be the value in which the 
primal solution is incremented in the i-th iteration. 



Let Dl(l) = J2e£E C eVe and D * (*) = T, 



rER • 



be the 



respective values of the first and second parts of the 
dual solution at the end of the i-ih iteration, and let 
D(i) = Di(i) + Z?2(i). Also, let D denote the value of 
the optimal solution for the dual linear program. 



• Let a(i) denote the normalized length of the path se- 
lected after the end of the i-th iteration. Note that 
if path p is selected in the (i + l)-th iteration then 

°W = ^bl = ^E eep 2/e- 
Correctness and Truthfulness. The following lemmas 
establish the feasibility of the solution, and the monotonicity 
and exactness of the algorithm. 

Lemma 3.3. Algorithm Bounded-UFP(e) outputs a feasi- 
ble solution. 

Proof. Assume by contradiction that the output of the 
algorithm is not feasible. Let p be the first path that induces 
a violation in the capacity of edge e in the ^-th iteration, and 
let V be the family of paths selected before the £-th iteration, 
which consist of e. Since p induces a capacity violation then 
E„g-p dp + dp > c e . Clearly, since dp G (0, 1] it follows that 
E e -p dp > c e — 1. Consequently, we get that 



t-i _ o 



> e 



n 






>e £ 



a e(B-l) 



where the last inequality results from the fact that ^— - is 
an increasing monotonic function for all x > 1, and since 
Ce > B > 1. Inspecting the main loop stopping condition, 
i.e. line 5 in the algorithm, it follows that the algorithm had 
to exit the loop. This implies that the algorithm could not 
have executed the £-th iteration and thus, could not have 
selected p, a contradiction. □ 

Lemma 3.4. Algorithm Bounded-UFP(e) is monotone and 
exact w.r.t. to the demand and value of every request. 

Proof. Consider a request r selected to be routed using 
path p r in the £-ih iteration of the algorithm, which has a 
respective demand and value of d r and v r . Now, suppose 
that r had a demand of d r < d r , a value of v r > v r , and 
the demands and values of all the other requests were fixed. 
For the sake of monotonicity, we need to prove that the al- 
gorithm would have selected r in the latter case, i.e. when 
its demand and value were d r and v r , respectively. If r is 
selected by the algorithm in the first £—1 iterations then we 
are done. Otherwise, lets consider the ^-th iteration. One 
can easily observe that in the first £ — 1 iterations of the 
algorithm, the same set of requests is selected to be routed 
using the same set of paths whether the demand and value 
of r is (d r ,v r ) or (d r ,v r ). Respectively, the same set of un- 
selected requests remain. Note that 4p < — . Hence, since 

the path p r minimizes -£ E e p 2/ e over an y P a ^h P, which 
corresponds to an unselected request, when the demand and 
value of r is (d r ,v r ), so it does when the demand and value 
of r is (d r ,v r ). This implies that r must be selected by the 
algorithm in the £-ih iteration. 

The exactness of the algorithm is clear, as the algorithm 
may route the exact demand of every request selected, and 
may not route anything otherwise. □ 

Approximation. We now turn to prove that the algorithm 
yields an approximation ratio that approaches -^ri- We be- 
gin by stating three claims, which will be utilized later. 

Claim 3.5. An increasing sequence {ao,a±,. .. ,a t +i,a} 



satisfies J2 i=0 



i+l— °H 



< In I 



Proof. For every < i < t, 

< In = ln(a - on) - mia - a i+1 ) , 

a — an \a — cti+i / 

where the inequality follows from ln(l + x) < x by substi- 
tutmg x = — — 



Accordingly, this implies that 



£( £j 5^)*£(^«-<*)- h *»-«*o) 



ln(a — qq) — In (a — Ott+i) — In 



a — ctQ 
a — ctt+i 



D 



Claim 3.6. a(i) < D _ ^ , \n , ^ n every iteration i > 0. 

Proof. Consider the (i + l)-th iteration. Let p denote 
the path that is selected in this iteration. The path p cor- 
responds to an unselected request such that -^ X^e? Ve is 
minimal. Namely, every other path p' , which corresponds 
to another unselected request, satisfies 



d 



d 



;J- 5>« ^ £ tf = ««. thus d " £ ^ - v ■ 

P egp' v e£p egp' V ' 

This implies that if we multiply yl by a(i)~ , for every 
e 6 E, then all the dual linear program constraints become 
satisfied. Consequently, the set of variables (j/'a(i) _1 ,2 l ) 
constitutes a feasible fractional solution to the dual linear 
program and therefore, D < D 1 (i)a(i)~ 1 + D2(i). □ 

Claim 3.7. Di(i+1) < Di(i)+Be(l+e)-AP(i+l)-a(i), 
for every i > 0. 

Proof. Consider the (i + l)-th iteration. Let p denote 
the path that is selected in this iteration, and let d p and 
v p denote its respective demand and value. Inspecting the 
algorithm, one can derive that 



£ c e yi +1 = £ Ce ^ + £ Ce ^ ■ e "" 

eG-B e€E eGp 

e<£p 

<E^+E^(i + ^ + (^) 2 ) 

eS-B egp \ e e / 

- £ Ce ^ + Be ( 1 + £ ) d P £ ^ 
eGB eGp 

= £ Ce ^ + B < 1 + £ ) ' AP ^ + X ) ' °(*) ■ 

eG-B 

The first inequality is due to the fact that e a < 1 + a + a 2 for 
any a £ [0, 1], and the fact that p £ (0, 1]. The second 
inequality holds since 



eBd v 



eBd v \ 2 



+ (-^ j j < ^Bd p + e 2 Bd v = Be(l + e)d v , 

(1) 

where the inequality in (1) follows from the observations 
that dp < dp, and — < 1. Finally, the last equality fol- 
lows from the definition of a(i), which can be rewritten as 
^p XeGp Ve = u p a (*)> an< i ^ ne observation that v p is the value 
in which the primal solution is incremented in the (i + l)-th 
iteration. Recalling that Di(i) = Xbgb; c e yl completes the 
proof. □ 



We are now ready to establish the approximation guarantee 
of the algorithm. 

Lemma 3.8. Algorithm Bounded-UFP(e) returns an ((i+ 
6e)^j) -approximate solution for the -Ap -bounded unsplit- 
table flow problem, for any e G (0, I]. 

Proof. One can easily notice, by inspecting the stoping 
condition of the main loop, that when the algorithm termi- 
nates, either £ = or Y.eeE c <=V<= > e e(s ~ 1) . If C = then 
it follows that the algorithm succeeded to satisfy all the re- 
quests and thus, its output is optimal. Consequently, in the 
remainder of the proof, we shall assume that X^gb: Ce V e ^ 
e eCS_1) . For every i > 0, 



Z>i(i + 1) < Dx(i) + Be(l + e) ■ AP(i + 1) ■ a(i) 



< D 1 (i)+Be{l + e)-AP(i + l) 



Di(i) 



D~D 2 {i) 



Di(t)(l + Be(l + e) 



AP(i + 1) 
D-D 2 (i) 



l( i )e (*«^mm: 



< D 

where the first and second inequalities follow from Claim 3.7 
and Claim 3.6, respectively, and the last inequality is due to 
the fact that 1 + a < e a . This implies that 

Di(i + 1) < Di(0)e( SE(1+e)E 5=° o-otu)) 

^ e (He= + Be(l +e)E j =0 ^^)) 

where the first inequality results from the expansion of the 
recursion, and the second one follows from -Di(O) < e Be , 
which is obtained by noticing that -Di(O) = Xegb; Ce Ve = m > 
and recalling that B > ^-p- by definition. Lets assume that 
the algorithm terminates after t+ 1 iterations. Accordingly, 
using our prior assumption that Xegb; CeVe > e"^ B ~ ', we 
get that 



a <B-l) 



<Di(*+l) <e 



( flB » + B«(i + e)Ej-o^m) 



One can validate that Dz(j) = P(j) by inspecting the vari- 
ables alterations in line 12 of the algorithm, and their af- 
fect on -D2O) and P(j)- Hence, we derive that D _jyr \ — 
jj_p, ■» ■ Consequently, we can apply Claim 3.5, while 
recalling that P(0) = and P(t + 1) = P, and yield 

e{B - 1) < Be 2 + Be(l + e) In (■£?) . 



Because i=2s < t(B -}}7, 



any positive e, we attain 1 — 3e < In 
simplified further to give ■£ < 



and since 1 — 3e < \ 2e for 
D ^ This can be 



D-P ) 

p _= ,,(1-30-! + ! ^ a+««)i=T- 

Recall that D is the value of the optimal solution for the 
dual linear program and thus, using the weak LP duality 
completes the proof. □ 

Proof of Theorem 3.1. The correctness of the algorithm 
is due to Lemma 3.3, the approximation guarantee is estab- 
lished in Lemma 3.8, and the monotonicity and exactness 
are presented in Lemma 3.4. Finally, it is clear that the 
running time of the algorithm is polynomial, and in fact, 
if we denote the number of requests by \1Z\, one can easily 
validate that the number of iterations is bounded by \1Z\, 
and every iteration takes time proportional to \R\ shortest 
path computations. ■ 



3.3 Inapproximability result 

In the following, we introduce two input instances for the 
problem under consideration that lower bound the perfor- 
mance guarantee of any algorithm, which is part of the rea- 
sonable iterative path minimizing algorithms family. Specif- 
ically, the first input instance proves that any such algo- 
rithm cannot yield an approximation guarantee better than 
— rj — o(l). This demonstrates that the analysis of algorithm 
Bounded-UFP is tight, and that it is the "best" algorithm in 
the aforementioned family of algorithms. The second input 
instance establishes a lower bound of | on the approxima- 
tion ratio of any such algorithm for this problem in its ut- 
most generality, i.e. when the underlying graph is undirected 
and the luinimal edge capacity is arbitrarily large. In partic- 
ular, this suggests that even if we ease the problem setting, 
e.g. assume that the minimal capacity of an edge is Q(m) 
instead of f2(lnm), no algorithm in the aforesaid family can 
achieve PTAS. 

Prior to describing the finer details of our approach, we 
introduce the notion of a reasonable function, which is a key 
ingredient in the definition of a reasonable iterative path 
minimizing algorithm. Note that reasonable functions have 
a similar flavor to the min functions introduced by Archer 
and Tardos [3]. Nevertheless, they are still quite different. 
Let S = Ure^iP : P is a simple path between s r and t r }. 

Definition 3.9. Let g : S — * R be a function, which as- 
signs an arbitrary priority to every path. Such a function 
is called reasonable if under the assumption that the capac- 
ities of all the edges are identical, and both the demand and 
value of every request are unit, it follows that g{q) < g{q ) 
for any valid unsplittable flow, and any two paths q, q' 6 S 
that satisfy 



q consists of k edges, q' consists of k! 
k'. 



and k < 



• fi < f[, for every 1 < i < k, where (/1, f 2 , . . . , f k ) 
and (/x, / 2 , . . . , fy) are non-increasing vectors, which 
indicate the flow routed through the edges of q and q 
with respect to the valid unsplittable flow. 

Definition 3.10. An algorithm is referred to as reason- 
able iterative path minimizing algorithm, if it iteratively se- 
lects a path that minimizes a reasonable function over all the 
paths that correspond to unselected requests. 

One can verify that algorithm Bounded-UFP minimizes 



the function h(p) — — J^ 



e£p c e 



where f e denotes the 
flow routed through edge e. This function is reasonable since 
when we assume that both the demand and value of every 
request are unit, and the capacities of all the edges are iden- 
tical, say B, then it reduces to h(p) = jj X^ep eE ^% which 
clearly satisfies h(q) < h(q'), for any two paths q and q' 
that meet the properties indicated in Definition 3.9. Con- 
sequently, algorithm Bounded-UFP is a reasonable iterative 
path minimizing algorithm. We note that reasonability cap- 
tures a broad class of functions. For example, the reasonable 
function h\(p) = ln(l + |p|) • h(p) is similar to the func- 
tion used by algorithm Bounded-UFP but it is mildly biased 
towards paths with less edges. Another example of a rea- 
sonable function is /12 (p) — — Yl e€ — although it is not 
clear why anyone would like to use it. We are now ready to 
establish the main result of this subsection. 



Theorem 3.11. The approximation ratio of any reason- 
able iterative path minimizing algorithm for the unsplittable 
flow problem when the graph is directed cannot be better than 
'^ T -o(l) ; /orB = (7n 1 '' 2 ). 




Figure 2: A directed graph in which every vertex 
Si has a directed edge to every vertex Vj such that 
j > i. Additionally, the capacities of all the edges 
are identical and equal to B. 

Proof. Suppose we are given the directed graph G — 
(V, E) schematically described in Figure 2, and the set of 
requests is 

ft = {(ai,t,l,l),(a3,t,l, !),...,(«,*, 1,1)} • 



B requests B requests 



B requests 



In order to simplify the presentation and analysis of the 
lower bound instance, we introduce the following assump- 
tion, which will be tackled later. We assume that when 
there is more than one path, which minimizes the value of 
the reasonable function used by the iterative path minimiz- 
ing algorithm, the algorithm selects one of them arbitrarily. 
Accordingly, we premise that it selects a path (si,Vj,t) in 
which i is minimal, and j is maximal with respect to all 
the minimizing paths that their source vertex is s;. For ex- 
ample, in the initial B iterations of the algorithm, all the 
requests that their terminal vertices are (si,t) are satisfied 
using paths that use vertices Vj such that j — £,... , l — B + 1. 
In the subsequent B iterations, all the requests that their 
terminal vertices are (s2,t) are satisfied using paths that use 
vertices Vj such that j = £ — B, . . . , £ — 2B + 1, and so on. 
Simulating the execution of a reasonable iterative path min- 
imizing algorithm, while ignoring integrality issues that will 
be resolved later, we get that 

• For any integer 1 < q < B, at the end of the first 
^X^r=i("BTi) r iterations, all and only the requests, 
whose terminal vertices are (si,t) such that i < I ■ 
(1 — (-gfrj") 9 )i are satisfied, and all the {vj,t) edges 
such that j > £ • (1 — (-grrj-) 9 ) have a flow load of q or 
equivalently, a residual capacity of B — q. 

• After ^5Z r=1 (-gTTf) r iterations, the algorithm cannot 
route more requests and thus, it stops 2 . 



2 An algorithm might stop even sooner, e.g. algorithm 
Bounded-UFP stops after the while condition fails. How- 
ever, analyzing the case that the algorithm stops when it 
cannot route more requests just affirms the lower bound. 



Consequently, since all and only the requests, whose ter- 
minal vertices are (st,t) such that i < I ■ (1 — (-gqrj-) 5 ), 
are satisfied when the algorithm stops, it follows that the 
value of the solution that the algorithm outputs is at most 
B£.(l-(^ T ) B ) = B£-(l-(l- 5 i T ) B )<B£.(l-I). On 
the other hand, an optimal solution clearly has a value of 
B£, e.g. route every request of the form (si,t, 1, 1) through 
the directed path (si,Vi,t). Thus, we get that the approxi- 
mation ratio of the algorithm cannot be better than ^§y . 

We now drop the integrality assumption. Namely, in the 
above analysis, we have assumed that ^X^r=i("BTr) r * s m t e_ 
gral, for any integer q. This clearly may not be true. How- 
ever, one can resolve this issue by applying a more care- 
ful analysis, and yield that only the requests, whose ter- 
minal vertices are (si,t) such that i < £ ■ (1 — (-gfrf) 9 ) + 

5^fe=o("BTr) fe > may become satisfied. Since 5^fc=o {Sp[) k < 
B, it follows that the value of the solution that the algorithm 
achieves might increase by no more than B 2 , in respect to 
Bl ■ (1 — -). Since the number of edges in the graph is 

m = t + X]fe=i & — 2^ 2 , and B = olm 1 ' 2 ), we obtain that 
B — o(£). Consequently, the analysis of the lower bound de- 
grades by at most o(l), i.e. the approximation ratio of any 
algorithm cannot be better than — §j- — o(l). 

Next, we tackle the decisions assumption. Specifically, we 
have assumed that the decision of any algorithm between 
same valued minimizing paths is arbitrary and hence, one 
may ask if a specific tie-breaking rule can lead to better 
results. We can resolve this matter by constructing a similar 
input instance, which forces any algorithm to make similar 
"bad" decisions. Essentially, one way to achieve it is to 
replace every (si, Vj) edge by a directed path with ii-\- 1 — j 
edges. The reason that an algorithm makes "bad" decisions 
on this instance dues to the reasonability property, i.e. any 
reasonable algorithm "prefers" paths with less edges. Note 
that this instance supports the same lower bound, but has 
a somewhat stricter constraint on the value of B, i.e. since 



= 0(£ ), B needs to satisfy B = o(? 



,V^ 



□ 



The next theorem demonstrates that even if we ease the 
problem setting, no reasonable iterative path minimizing al- 
gorithm can achieve PTAS. 

Theorem 3.12. The approximation ratio of any reason- 
able iterative path minimizing algorithm for the unsplittable 
flow problem cannot be better than ^, for any B, and even 
when the graph is undirected. 




Figure 3: An undirected graph in which the capaci- 
ties of all the edges are identically equal to B. 

Proof. Suppose we are given the undirected graph G — 
(V,E) schematically described in Figure 3, and the set of 
requests is 

71 = { {vi,V 3 , 1, l) ; frt,tt6, 1, 1) ; frl,^ 6 , 1, l) ; j>3,tt4, 1, 1) } • 
B requests B requests B requests B requests 



Clearly, an optimal solution for this instance has a value of 
4B, e.g. route every request of the form (vi, Vz, 1, 1) through 
the path (t>i, ^2,^3), any request of the form («4, Vq,1, 1) 
through the path («4, V5,ve), all the requests of the form 
(vi,ve, 1, 1) through the path (vi,V7,Ve), and every request 
of the form (113,114,1,1) through the path (1)3,1)7,1)4). 

Simulating the execution of a reasonable iterative path 
minimizing algorithm in the initial four iterations, one can 
easily validate that the algorithm may select the four paths 
(vi,vr,V3),(v4,V7,Ve), («i,V2,«3), and (v4,vs,v 6 ) since each 
one of these paths is one of the minimizing paths in the 
corresponding iteration, and by that satisfy two (vi, V3, 1, 1) 
requests and two (v4,vq, 1,1) requests. In addition, notice 
that at the end of this four iterations phase, every edge has 
a residual capacity of B — 1. Arguments similar to those 
used in this initial four iterations phase can be applied in 
another -j — 1 phases, each of four iterations, to demonstrate 
that the algorithm acts exactly the same. Consequently, 
after -|- phases, every edge has a residual capacity of -j, all 
the (1)1,1)3, 1, 1) requests and the (v4,ve, 1, 1) requests were 
satisfied, and the remaining requests are 

n = {(v 1 ,V 6 ,l,l),(v 3 ,V 4 ,l,l)} . 



B requests 



B requests 



In this current state, any algorithm can satisfy at most B 
requests from TZ. This is the result of the fact that any path 
from vi to vq and any path from V3 to V4 must use either edge 
(1)1,1)7) or edge («3,«7), and the fact that the total residual 
capacity of these edges sums to B. Thus, the solution that 
the algorithm outputs has value of at most SB. □ 

Corollary 3.13. No reasonable iterative path minimiz- 
ing algorithm for the unsplittable flow problem can yield a 
PTAS. 

4. SINGLE-MINDED MULTI-UNIT 
COMBINATORIAL AUCTION 

4.1 The algorithm 

In this subsection, we design a deterministic monotone 
algorithm for the f2(lnm)-bounded multi-unit combinato- 
rial auction problem, whose approximation ratio approaches 
— §j. We begin by demonstrating that the single-minded 
multi-unit combinatorial auction problem can be formulated 
as a simplified special case of the integer linear program of 
the unsplittable flow problem, and then we turn to special- 
ize algorithm Bounded-UFP for the problem under consider- 
ation. 

The single-minded multi-unit combinatorial auction prob- 
lem can be formulated as a special case of the integer linear 
program of the unsplittable flow problem by letting S r to 
denote the singleton set of U r , i.e. S r = {U r }, and replacing 
e, E and d„ in the integer linear program of the unsplit- 
table flow problem with u, U and 1, respectively. Similarly, 
the relaxation of the integer linear program of the single- 
minded multi-unit combinatorial auction problem, and its 
dual can also be formulated as special cases of the corre- 
sponding linear programs. Consequently, the primal-dual 
algorithm Bounded- MUCA, formally described below, is a 
specialized version of algorithm Bounded-UFP, in which the 
path selection procedure, i.e. lines 6-8 in algorithm Bounded- 
UFP, was neglected, and the demand terms were omitted. 



Algorithm 2 Bounded-MUCA(e) 



Input: An accuracy parameter e £ (0, 1] 
Output: A set W of requests to be satisfied 

1: Let C be a list of all the requests, and let W be an empty 

set 
2: for all u 6 U do y u — — end for 

3: while (C / and £ uGC/ c uJ/u < e^- 1 )) do 
4: Let f be the request, which minimizes — X^ugu V u 
with respect to every r £ C 



for all u £ Uf do y u ~ y u ■ e 



B/c 



end for 



Add f to W, and remove f from £ 
end while 

return W 



Theorem 4.1. T/ie algorithm Bounded-MUCA(^) returns 
a feasible {(i+n)-^-^)- approximate solution for the fl(^-p-)- 
bounded multi-unit combinatorial auction problem, for any 
e 6 (0,1], runs in polynomial-time, and is monotone and 
exact w.r.t. the value of each request. 

Proof. Since algorithm Bounded-MUCA is a simplified 
version of algorithm Bounded-UFP, the analysis of Theo- 
rem 3.1 also applies in this case. □ 

It is worth noting that the algorithm can even be em- 
ployed to a generalized version of the problem in which the 
demand of every request, i.e. the desired bundle of items, 
is part of the type of the request and therefore, agents may 
be dishonest about it. Note that in this case, monotonic- 
ity can be easily established by arguments similar to those 
used in Theorem 3.4, and the additional observation that 
J2ueu r J/« ^ H2ueu t Vu> for an y Ur Q U r . Remark that this 
setting is referred to as the unknown single-minded case [14]. 

Corollary 4.2. For all values e e (0,1], there exists 
a polynomial-time truthful ((i+e)-^j) -approximation mecha- 
nism for the f2(ln m) -bounded multi-unit combinatorial auc- 
tion problem among unknown single-minded agents. 

4.2 Inapproximability result 

We demonstrate that any reasonable iterative bundle min- 
imizing algorithm cannot achieve an approximation ratio 
better than |. Essentially, this suggests that no algorithm 
in the aforesaid family can achieve PTAS. 

We start by formally defining the notion of a reasonable 
iterative bundle minimizing algorithm. We remark that the 
following definitions are just a refinement of the definitions 
made for the unsplittable flow problem. Let S = {J reK {U r }. 

Definition 4.3. Let g : S — * R be a function, which as- 
signs an arbitrary priority to every bundle. Such a function 
is called reasonable if under the assumption that the multi- 
plicities of all the items are identical, and the value of every 
request is unit, it follows that g(T) < g(T') for any valid 
multi-unit allocation 3 , and any two bundles T,T' € S that 
satisfy 



3 A valid multi-unit allocation can be succinctly described 
as an allocation of non-identical items to requests such that 
the number of allocated copies of every item does not exceed 
its multiplicity. 



• T consists of k items, T consists of k items, and k < 
k'. 

• .A < ft, for every 1 < i < k, where (/i,/ 2 , . . . ,/*) 
and {f i, f %,■■■, f'y) are non-increasing vectors, which 
indicate the number of allocated copies of the items of 
T and T with respect to the valid multi-unit allocation. 

Definition 4.4. An algorithm is referred to as reason- 
able iterative bundle minimizing algorithm, if it iteratively 
selects a bundle that minimizes a reasonable function over 
all the bundles that correspond to unselected requests. 

Note that algorithm Bounded-MUCA minimizes the func- 
tion h(s) = — X^uss — e c " ' wnere f u denotes the number 
of allocated copies of item u. One can easily argue that this 
function is reasonable and thus, algorithm Bounded-MUCA 
is a reasonable iterative bundle minimizing algorithm. We 
are now ready to prove the core result of this subsection. 

Theorem 4.5. The approximation ratio of any reason- 
able iterative bundle minimizing algorithm for the single- 
minded multi-unit combinatorial auction problem cannot be 
better than |. 

Proof. Let m be a multiple of p ■ (p + 1), where p > 
3 is a constant odd integer. Suppose we are given a set 
U of m items such that the multiplicities of all the items 
are identical and equal to B, and let Uf =1 U?=i Ui,j be a 
partition of U into p ■ (p + 1) disjoint sets, each consists of 
,r m +1 -. items. Additionally, suppose that 1Z consists of unit 
value requests of two types: 

D_ 
2 

for every 



1. -j requests that consist of the items Ut — U?=i Ut,h 



2. -j requests that consist of the items Ui,2t-i U Ui,2t U 
U?=2 Ui,2t-i, for every £ — 1, . . . , £±-, and -| requests 
that consist of the items Ui^e-i U Ui,2i U {J^ =2 ^i,2ti 



for every I = 1, 



p+1 
2 ■ 



Figure 4 schematically describes a concrete input instance. 

One can verify that the an optimal solution for this in- 
stance has a value of pB, for example by selecting all the 
requests except for the §■ requests that consist of Ui. 

Simulating the execution of a reasonable iterative bun- 
dle minimizing algorithm, one can easily validate that the 
algorithm may incrementally select all the requests of the 
first type, e.g. it may repeatedly select a request that con- 
sists of Ui, then a request that consists of U2 and so on 
until Up. Consequently, after -j phases, each of p steps, 
all the requests of the first type are satisfied, the value of 
the current solution is \B, and every item has a residual 
multiplicity of -5-. Notice that U\ consists of — items, and 



any request of the second type consists of 



items of 



p-(p+i) 

U\. Hence, by simple counting arguments it follows that in 
the current state, any algorithm cannot satisfy more than 
(f ' ir)/ P -( 2 p+i) = E i LB requests of the second type. There- 
fore, we obtain that the value of the solution that any rea- 
sonable iterative bundle minimizing algorithm outputs is no 
more than p ^ B and thus, as p tends to infinity, the inap- 



proximability ratio approaches 



□ 



Corollary 4.6. No reasonable iterative bundle minimiz- 
ing algorithm for the single-minded multi-unit combinatorial 
auction problem yields a PTAS. 
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Algorithm 3 Bounded-UFP-Repeat(e) 



1 



Figure 4: The input instance in the case of p = 3, 
and m = 12. Note that the set of items Ui=i U,=i Ui,j 
are represented by the black dots on the left of the 
figure, and every collection of y requests is repre- 
sented by a vertical group of rectangles such that 
each rectangle consists of the set of items, whose 
representing dots are in the left-projection zone of 
the rectangle. Also note that the first type of re- 
quests are located left of the dotted line, whereas 
the second type of requests appear on its right. 



5. UNSPLITTABLE FLOW WITH 
REPETITIONS PROBLEM 

In this section, we study the f2(lnm)-bounded unsplit- 
table flow with repetitions problem. This problem is a vari- 
ant of the corresponding unsplittable flow problem in which 
we are allowed to satisfy every request multiple times using 
possibly multiple paths, and the profit gained is propor- 
tional to the number of times that every request is satisfied. 
In sharp contrast with our prior results, we demonstrate 
that this version admits a deterministic primal-dual based 
algorithm, which yields an (1 + ^-approximation. 

5.1 The algorithm 

In the following, wc contrive an (1 + e)-approximation 
algorithm, named Bounded-UFP-Repeat, for the fi(ln re- 
bounded unsplittable flow with repetitions problem, which 
is based on a primal-dual approach. The following theorem 
digests the properties of the algorithm. 

Theorem 5.1. Algorithm Bounded-UFP-Repeat(^) is an 
(1 + e) -approximation for the Q^^p-) -bounded unsplittable 
flow with repetitions problem, for any e G (0, 1], with running 
time polynomial in m and ""',''') . 

In Figure 5, we present the primal-dual formulation of the 
underlying problem. 

We now turn to analyze algorithm Bounded-UFP-Repeat. 
For the sake of simplicity, we shall use the same notation, 
which was introduced in Subsection 3.1, with one exception. 
Namely, since there are no z-type variables in the dual linear 
program, we let D(i) = Xleefi c eVe denote the value of the 
dual solution at the end of the i-th iteration, and neglect 
D\(i) and D2(i). Notice that Lemma 3.3 is applicable also in 
this case and thus, the correctness of the algorithm follows. 
Consequently, in the sequel, we prove that the algorithm 
under consideration achieves an approximation ratio of (1 + 



Input: An accuracy parameter e G (0, 1] 
Output: A (request, path) pairs multiset W, which 
holds the requests to be allocated 

Let £ be a list of all the requests, and let W be an empty 

multiset 

for all e G E do jy e = — end for 

< e^-V) 



while (Z eeB c e y e <e^- L >) do 
for every r € C do 

Let p r be the shortest path between s r and t r in 
G with respect to the weights y e , and 
let \p r \ = Eeg Pr y<= be its length 
end for 

Let f be the request, which minimizes — \p r 
respect to every r € £ 



with 



8 

9 

10 


for all e£jif do y 
Add (r,Pf) to W 
end while 


11 


return VV 



Ve 



eBdf/c 



end for 



max 


E w -- 

r6R 


BES r 


■) 


s.t. 


E * 

s£S\e£s 


sCls J^ Ce 


VeG-E 




i,6N 




VsG S 


min 


E c< * 

eS-E 






s.t. 


dr E ^ e " 


IV W G !Z,\/s G S r 




e£s 








y e >0 




VeG E 



Figure 5: The integer linear program of the unsplit- 
table flow with repetitions problem (top), and the 
dual of its relaxation (bottom). Note that S r denotes 
the set of all the simple paths between s, and t r in 
G, S = {Jr^-^Sr, and d s and v 3 denote the respective 
demand and value of path s. 



e). We begin by establishing a analogous claim to Claim 3.6, 
which upper bounds a(i). 

Claim 5.2. a(i) < —$■, in every iteration i > 0. 

Proof. Consider the (i + l)-th iteration, and let p denote 
the path, which is selected in this iteration. The path p 
corresponds to a request such that -^- ^2 e& yl is minimal. 
Namely, every other path p' , which corresponds to a request, 
satisfies 



dpi 
v 



d 



— Y] We > — Y] 2/e = ot(i), thus d p , y] 



He 



eGp' 



eGp 



e£p 



', ot(i) 



> v„ 



This implies that if we multiply y' e by a(i) _1 , for every 
e € E, then all the dual linear program constraints be- 
come satisfied, i.e. the modified variables constitute a fea- 
sible fractional solution to the dual linear program. Hence, 
D < D(i)a(i) _1 . □ 



Next, we prove that the algorithm achieves the claimed 
(1 + e)-approximation for the problem under consideration. 

Lemma 5.3. Bounded-UFP-Repeat(e) returns an (l + 6e)- 
approximate solution for the ^^ -bounded unsplittable flow 
with repetitions problem, for any e £ (0, g]. 

Proof. For every i > 0, one can derive that 

D(i + 1) < D(i) + Be(l + e) ■ AP(i + 1) ■ a(i) 

D(i) 



< D(i) + Bt{l + e) ■ AP(i + 1) 



D 



D(i)[l + Be{l + e 



AP(i + l) 



D 



< D 



(i)eM^) 



The first inequality follows from Claim 3.7, while noticing 
that D(i) may replace Di(i). The second inequality is due 
to Claim 5.2. Finally, the last inequality results from the 
fact that 1 + a < e a . This implies that 

Di(i + 1) < d^e^'^^ 11 ) 
< ( fl ^+B.(i+«)i:j =0 ^±ii) 



where the first inequality results from the expansion of the 
recursion, and the second inequality follows from -Di(O) < 
e Se , which dues to -Di(O) = m, and B > i^^. Lets assume 
that the algorithm terminates after t+1 iterations. Accord- 
ingly, inspecting the stoping condition of the main loop, we 
attain that 



<B-i) 



< Dl (t + 1) < e ( B < 2+B < 1+ ^U ^&ll) 



This can be simplified to 



e(B-V)-Be' 
Be(l+e) 



<£.U 



AP(j + l) 



No- 



tice that J]] ._ AP(j + 1) is a telescopic sum that is equal to 
P(i+1)-P(0) = P, where the equality dues to P(t+1) = P 
and -P(O) = 0. In addition, notice that I — 3e < WTiT^T 
for any positive e. Consequently, we obtain 1 — 3e < -£. Sim- 
plifying this expression even further yields -^ < * < I + 
6e, which by the weak LP duality establishes the lemma. □ 

We are now ready to prove the main theorem of this sub- 
section. 

Proof of Theorem 5.1. As noted before, the correctness 
of the algorithm directly follows from Lemma 3.3. In addi- 
tion, the approximation guarantee is proved in Lemma 5.3. 
We now turn to argue that the running time of the algo- 
rithm is polynomial in m and f 3 ^, where c ma x = max e {c e }, 
and d m in = min r {d r }. Consider some edge e 6 E. Recall 
that J/e = — and y\ < — e eS , where t denotes the index of 
the last iteration of the algorithm. In addition, notice that 
every time that the algorithm increments y e , it is by a mul- 
tiplicative factor of at least e eSdmin ' Cmax . Consequently, the 
number of iterations in which y e is incremented is at most 
^ max . This implies that the running time of the algorithm 
is bounded by m ^ max , as there are m edges. ■ 
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